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ABSTRACT 



A comparison of a simple recursive ray acoustics algorithm versus a ray acoustics 
algorithm based on solving a system of first-order ordinary differential equations was 
conducted. The recursive ray acoustics (RRA) algorithm was found to be accurate and 
relatively fast. The RIC\ algorithm is capable of handling sound speed as a function of 
all three spatial coordinates, and this capability was demonstrated. 

Two separate methods of representing a sound-speed profile (SSP) based on data 
points were examined: Akima cubic spline and spatial Fourier series (SFS). The SFS 
representation encountered difficulties in accurately modeling SSPs. Various techniques 
were applied to improve the SFS sound-speed representation. While accurate sound- 
speed fits were eventually achieved, difficulties remained in the SFS modeling of first and 
second-order derivatives of the sound-speed data. 

The RRA algorithm was tested using the SFS sound-speed representation and found 
to be significantly inaccurate. 

A demonstration was conducted of the ability of the SFS sound-speed represen- 
tation to incorporate randomness in the SSP. 
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I. INTRODUCTION 



The primary focus of this thesis was to test the accuracy and capabilities of a re- 
cursive ray acoustics (RRA) algorithm that can be used to calculate the position, angles 
of propagation, travel time, and path length along a ray path and to generate ray trace 
plots for speeds of sound that are functions of all three spatial coordinates. The RICA, 
algorithm is an extension and slight modification of an algorithm briefly discussed by 
Klein [Ref. 1: pp. 29-31]. The algorithm discussed by Klein [Ref. 1: pp. 29-31 ] was e.x- 
pressed in terms of an index of refraction while herein, the RRA algorithm is expressed 
directly in terms of the speed of sound. No computer simulation results were presented 
by Klein [Ref I: pp. 29-31 ]. 

The accuracy of the RRA algorithm was tested by comparing its computer simu- 
lation results with those obtained from a previously developed ray acoustics algorithm 
based on solving a system of first-order ordinary difierential equations (ODE) [Ref 2j. 
The ODE algorithm as applied to ray acoustics was discussed in detail by Lim [Ref 3j. 
Both the ODE algorithm and the RRA algorithm used the IMSL Version 10 double 
precision computer program DCSAKM to obtain an Akima cubic spline fit to one- 
dimensional, depth-dependent, sound-speed data for testing and comparison purposes. 
The computer simulation results are presented in Chapter III, Section A. An Akima 
cubic spline was chosen since it is designed to produce a curve that matches the shape 
of the data while minimizing oscillations. In addition, the RR/\. algorithm was tested 
and shown to be capable of handling a sound-speed profile (SSP) that is a function of 
all three spatial coordinates. This capability was demonstrated by first separating the 
ocean into zones along the z axis. In each zone, the SSP is a different function, in gen- 
eral, of both cross-range x and depth y. A SSP with a gradient in the cross-range di- 
rection will cause a sound ray to leave its initial plane of propagation. Chapter III, 
Section B, includes demonstrations of the special three-dimensional capabilities of the 
RRA algorithm. 

A secondary objective of this thesis was to examine the use of a spatial Fourier series 
(SFS) to fit sound-speed data and to assess its utility relative to an Akima cubic spline 
fit. The appeal of the SFS lay in its representation of spectral content in terms of spatial 
frequencies, as well as its ability to represent derivatives. The precision of the two 
sound-speed fits, as well as that of the first and second-order derivatives of the sound- 
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speed data were compared since the first-order derivative as well as the sound-speed fit 
itself are utilized in the ray trace algorithms mentioned above. Chapter IV, Section A, 
contains the results of this comparison as generated by computer simulation. The SFS 
representation of a SSP was then employed within the RICA algorithm and tested by 
using the same SSPs that were used to compare the ODE and RIC\ algorithms. Chapter 
IV, Section B, discusses the accuracy of the RICA algorithm incorporating the SFS rep- 
resentation of the various SSPs. 

Additionally, random fluctuations in the SSP were modeled using the inherent suit- 
ability of a SFS representation for the inclusion of components at particular harmonics. 
This random element was intended to simulate fluctuations generated by such processes 
as internal waves or by inaccurate collection of sound-speed data. 
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II. THEORETICAL ANALYSIS 



A. THE RRA ALGORITHM 

If a sound source is located at Zq) and if a ray is launched with initial an- 

gles of propagation /?„ and 4>a (see Figure I ), then an estimate of the position of subse- 
quent points along the ray path is given by (see Figure 2 ) 

r/ = r/_i + Ar/.p /= 1,2,3,... (2.1) 

where 

Ar, = As,«,, (2.2) 

As, is an increment of arc length, the magnitude of which is at the user's discretion, 

ni = UiX + Viy + WiZ (2.3) 

is the unit vector along the ray path where 

Uj = sin /?,• cos 4>i , (2.4) 



V,- = cos /?,• , 



(2.5) 



and 



Wi = sin Pi sin <pi 



( 2 . 6 ) 



are the dimensionless direction cosines with respect to the x,y and z axes, respectively, 
and p, and 0, are the angles of propagation. 

Since, in general, is not identical with from Klein [Ref 1: pp. 29-31] we see 
that, as amplified and applied to the ocean acoustics environment in Ziomek [Ref 4: pp. 
225 and 230], 



t/(w(r)n(r)) 

ds 



= V/t(r), 



(2.7) 
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y 



Figure 1, Initial Angles of Propagation and (f>^ . 




Figure 2. Ray Geometry. 
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or, approximately, 



= V«(£;)A5, 



( 2 . 8 ) 



where «(£,) = h , , 



«(r/) = 



c{tri) 

c{ld 



(2.9) 



is the dimensionless index of refraction, c(£o) is the three-dimensional speed of sound at 

£o = +>'oi' + . (2-10) 



and, c(r,) is the three-dimensional speed of sound at 



A A A 

Ei = XiX +yiy +ZiZ , 



( 2 . 1 1 ) 



which is the position vector to a point on the ray path after the ith iteration. 
It follows from Equation (2.8) that 



and, as a result. 



w(r,) fii - rt(r,._j) «;_j = Vn{r,) As,- 



^ «(r, •_,)«/_, , ^Sl ^ ^ ^ 



( 2 . 12 ) 



(2.13) 



where, for better accuracy in the evaluation of the gradient, 



r/_, =r/_i +■ 



Asi 



7-1 A 






(2.14) 



which is the three-dimensional spatial location half-way along the arc length between 
r,-i and £„ and 



T7 ^ ^ , d A d A 



(2.15) 



is the gradient expressed in rectangular coordinates. It should be noted that although 
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A5, has been shown as a function of iteration number, it will be constant, in general, 
except when it needs to be adjusted to properly handle surface and bottom reflections 
and maximum down range limits. 

In terms of the direction cosines u, v, and w, 



and 



«(£/-! )?^/-l 


- + 


As 


^«(r*_i) 


(2.16) 


«(r/) 


«(r/) 




rt(r/_,)v;_, 




As 


^«(r/_i) 


(2.17) 


«(r/) 


- + • 


«(r,-) 




«(t:/-i)v^;- 


1 


As 




(2.18) 


«(r,) 


- + 


«(a) 


dz 



If c{r) = c{x,,y„z) is the three-dimensional speed of sound, and if Cq = c(£o) is the 
speed of sound at the source location, then it follows that 



Co c(r,) c(r,) 



n{ri) c(r;_,) Cq c(r,-_,) 



(2.19) 



and 



dn(r) 



n{r^ dx 



= - c(ri) 



dc(r) 






dx 



( 2 . 20 ) 



Analogous reasoning applied to the y and z coordinates yields results similar to 
Equation (2.20) , and substitution of these results into Equation (2.13) yields 



A C(i;) \ A 

^ c(/:,_i) 



I- c^(r) 



Vc(r) 



c = O-l 



( 2 . 21 ) 



6 



In addition, although Klein did not discuss travel time t and path length s calcu- 
lations, these quantities were computed using the following formulae; 



T = 




( 2 . 22 ) 



and 



s = ^A5,. (2.23) 

i=i 

Equations (2.1) through (2.23) constitute what is referred to herein as the recursive 
ray acoustics (RRA) algorithm first mentioned in Chapter 1. Input to this algorithm in 
the form of mathematical formulae for the sound-speed profile or some particular nu- 
merical fit to sound-speed data must be provided to enable the algorithm to function. 

B. SPATIAL FOURIER SERIES REPRESENTATION OF SOUND-SPEED DATA 

Given an ocean D meters deep with a one-dimensional, depth-dependent sound- 
speed profile (SSP), the speed of sound c(y) can be represented by the following spatial 
Fourier series (SFS): 



s 

cO) = ko I + 2^ I I cos{2nfy^y - Lc,), 0<y<D (2.24) 

5=1 



where |c,l and Lc, are the magnitude and phase(angle), respectively, of the complex 
Fourier series coefficient c, at harmonic s, kol is the magnitude of c, at harmonic s = 0 
and is the "DC component" or average value of the sound-speed data. 




(2.25) 



is the fundamental spatial frequency in the direction in cycles per meter where N is the 
total number of sound-speed data points taken and is the uniform sampling period 
in the y direction in meters per sample, and 5 is the highest harmonic. Derivatives of 
this Fourier series representation can be directly obtained. 
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The complex Fourier series coefTicients c, are obtained by taking the spatial discrete 
Fourier transform (DFT) of the sound-speed data as follows: 



v-i 



V 



IV 



■v< 



s = — 



n=0 



(2.26) 



where 



fVy = exp( +J2tzIN} 



(2.27) 



and, when sampling at the Nyquist rate, 

/V=25-M. (2.28) 



An attempt to represent a function by a truncated Fourier scries will always result 
in oscillations about the functions' true value. These oscillations normally have small 
enough amplitudes that they are of little consequence. However, an attempt to use a 
SFS representation of a linear SSP with a single constant gradient results in significant 
discrepancies in amplitude near the ocean surface and bottom, due to the inherent dis- 
continuity in the SSP that exists at these boundaries. This difficulty, known as "Gibbs' 
phenomenon", does not diminish as more harmonics, that is, sound-speed data points, 
are added to the Fourier scries. 

C. Lanczos [Ref 5] suggested that this phenomenon would be ameliorated if the 
Fourier series could be forced to converge more rapidly by modifying the Fourier series 
coefficients c, by a correction factor, such that the Lanczos smoothed Fourier series co- 
efi'icients are given by 



sin(s7r/S) 

snjS 



5 #= 0 , 



(2.29) 



where c, is given by Equation (2.25). Hamming (Ref 6] states that S in Equation (2.29) 
is usually replaced with S -I- I and, thus, our representation of the correction factor will 
be 



sm{snl{S + 1)) 
snIiS -M) ^ 



(2.30) 
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C INCORPORATION OF RANDOMNESS INTO THE SOUND-SPEED 
PROFILE 

The speed of sound as a function of depth, c(y) , can be represented as the sum of 

J 

deterministic and random components as follows: 

c(y) = coiy) + cn(y), (2.31) 

with Ci){y) being the deterministic and the random speeds of sound as functions of 
depth y. 

The autocorrelation function of Cj^(y) is , for any particulars^ andy , 

Rcp^y') = £^h0‘)40')l . (2.32) 

where c’^(y') is the complex conjugate of at depth y, 0 ^y ^ D . 

If 



C/tO) = exp( - j2nnfyj), 

n=-N 

for some set of Fourier series coefficients c„, and 



(2.33) 



(2.34) 



for D the ocean depth in meters, then 

c„ = -^ I c,(y) exp( +J2Knfy^y)ciy. (2.35) 

If Cfl(>’) is a random function of depth, then c, will be a random variable. 

Since 



ElcniCn-)*} = f f £{c^O)c^O')} 6xp( + j2nfyj^ny - n'y'))dydy' , (2.36) 

D *'0 *’0 



where 



E{cR(y)cR(y’)] = R,p,y’), 



(2.37) 
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and (q)‘ is the complex conjugate of c„. , and if it is assumed that CrO/) is a wide-sense 
stationary ( i.e., statistically homogeneous ) random process, then 

Rcp^y') = ^Ay), (2.38) 

where Ay =y — y' . If is represented by a finite Fourier series 

,v 

^CR(Ay) = Yj . (2.39) 

m=— /V 

then, as developed in Thomas [Ref. 7: pp. 148-153], 

~ 2iV -I- 1 (2.40) 

where is the power spectral density of Cs(y). 

Thomas [Ref. 7: pp 148-153] goes on to show that 

(2.41) 

where (3„„. is the Kronecker delta function. Therefore, 

I I } — ’'n— 2N + 1 (2.42) 

Hence, if c„ is a zero mean Gaussian random variable with variance r„ , that is, if 

c, = N{0,r„), (2.43) 

then Cr(v) will be a zero mean Gaussian random process. 

Therefore, 



\(y< y ) = ^c«(Ay) = ®^P( “ J'^^fya^y)^ 

n 

which represents the Inverse Discrete Fourier Transform of the power spectral density. 

Thus, knowledge of the power spectral density of Cfi(y) enables calculation of r„ via 
Equation (2.42), which in turn yields c„ via (2.43) and c^ly) via (2.33). 
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III. COMPUTER SIMULATION RESULTS 



A. RAY ACOUSTICS - COMPARATIVE ANALYSIS OF ODE VERSUS RRA 
ALGORITHMS 

Four distinct sound-speed profiles (SSPs) were examined for purposes of compar- 
ative analysis. For all cases, a sound source spatial location of Xq = 0 m and = 0 m 
was used. The source depth was established at a point half-way between the ocean 
surface and bottom, 100 m in most cases. A constant ocean depth of 200 m was mod- 
eled, unless otherwise noted, and for all cases three separate rays were examined. All 
rays were propagated initially in the direction of the positive z a.xis, that is, = 90°, but 
^0 was altered in order to generate unique rays. For SSPs 1 to 3, the launch angle /?„ 
was set to 45°, 85°, and 135° . SSP 4 was designed to simulate a SOFAR channel. 
Launch angles were therefore chosen to result in the first two rays propagating down 
range without reflection, i.e., /?„ = 85° and 95°, while the third launch angle was set to 
/Jo = 135° to test reflection as well. 

The four SSPs considered are defined as follows; 

• SSP 1: Constant speed of sound ( c = 1500 m, s ) 

• SSP 2: Linear SSP with a positive gradient ( c(y) = c(0) -I- gy m/s where the speed 
of sound at the surface c(0) = cO/ = 0) = 1500 m/s and g = -t-0.017 1/s ) 

• SSP 3: Linear SSP with a negative gradient ( c(y) = c(0) -t- gy m/s where the 
speed of sound at the surface c(0) = c(y = 0) = 1500 m/s and g = —0.017 1/s ) 

• SSP 4; Parabolic SSP with vertex at j;o where c(vo) = 1490 m/s and the speed of 
sound at both the ocean surface and bottom is equal to 1500 m/s 

( c(v) = c(yo) + a(y — >’o)^ ni/s with a = 0.001 ocean depth of 200 

m). 

The above four SSPs generated nine sound-speed data points which were input to 
the Akima cubic spline routine. The results from the ODE algorithm were used 
throughout as the benchmark for comparison purposes. The ODE algorithm used the 
IMSL Version 10 double precision computer program DIVPBS in order to solve the 
system of first-order ordinary differential equations. A range step size of 2.0 m and an 
error tolerance of 1 x 10~^ was used in DIVPBS, except for SSP 3, Ray 1 in Table 3, 
where the tolerance had to be reduced to 1 x 10”‘ in order for the ODE algorithm to run. 
The arc length step size used in the RRA algorithm was 2.0 m unless otherwise noted. 
The results of the various simulation runs are presented in the following seven tables. 
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Table 1. RAY ACOUSTICS RESULTS FOR SSP 1, RAYS 1 THROUGH 



Method 


J'o (m) 


fio (<lcg) 


r (km) 


y (m) 


/? (deg) 


T (sec) 


5 (km) 


ODE 


100.0 


45.0 


10.0 


100.00 


45.000 


9.428090 


14.1421 


RRA 


100.0 


45.0 


10.0 


99.95 


45.000 


9.428091 


14.1421 


ODE 


100.0 


85.0 


10.0 


174.89 


85.000 


6.692132 


10.0382 


RR^\ 


100.0 


85.0 


10.0 


174.88 


85.000 


6.692132 


10.0382 


ODE 


100.0 


135.0 


10.0 


100.00 


135.000 


9.428090 


14.1421 


RRA 


100.0 


135.0 


10.0 


100.05 


135.000 


9.428091 


14.1421 



For SSP 1 (constant speed of sound), Table 1 shows that at a horizontal range 
r = 10 km, the differences between the two algorithms are as follows: 

• between 0.01 and 0.05 m in depth >' and 

• 1 fi sec in travel time t. 
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Table 2. RAY ACOUSTICS RESULTS FOR SSP 2, RAYS 1 THROUGH 3 



Method 


>0 (ni) 


/?0 (deg) 


r (km) 


y (m) 


(deg) 


T (sec) 


s (km) 


ODE 


100.0 


45.0 


10.0 


99.99 


45.000 


9.417413 


14.1421 


RRA 


100.0 


45.0 


10.0 


99.95 


45.000 


9.417422 


14.1421 


ODE 


lOO.O 


85.0 


10.0 


162.50 


85.487 


6.683492 


10.0374 


RRA 


100.0 


85.0 


10.0 


162.57 


85.487 


6.683496 


10.0374 


ODE 


100.0 


135.0 


10.0 


100.01 


135.000 


9.417413 


14.1421 


RRA 


100.0 


135.0 


10.0 


100.05 


135.000 


9.417422 


14.1421 



For SSP 2 (linear SSP with a positive gradient), Table 2 shows that at a horizontal 
range r = 10 km, the differences betw’een the two algorithms are as follows: 

• between 0.04 and 0.07 m in depth j; and 

• between 4 and 9 fj. sec in travel time t. 
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Table 3. RAY ACOUSTICS RESULTS FOR SSP 3, RAYS I THROUGH 3 



.Method 


3o (m) 


Po (deg) 


r (km) 


y (m) 


/? (deg) 


T (sec) 


s (km) 


ODE 


100.0 


45.0 


10.0 


99.99 


45.000 


9.438784 


14.1421 


RRA 


100.0 


45.0 


10.0 


99.95 


45.000 


9.438795 


14.1421 


ODE 


100.0 


85.0 


10.0 


168.89 


84.512 


6.699404 


10.0380 


RRz\ 


100.0 


85.0 


10.0 


168.96 


84.511 


6.699409 


10.0380 


ODE 


100.0 


135.0 


10.0 


100.01 


135.000 


9.438784 


14.1421 


RRz\ 


100.0 


135.0 


10.0 


100.05 


135.000 


9.438794 


14.1421 



For SSP 3 (linear SSP with a negative gradient), Table 3 shows that at a horizontal 
range r = 10 km, the difTerences between the two algorithms are as follows: 

• between 0.04 and 0.07 m in depth y , 

• 0.001 deg in angle of arrival /? , and 

• between 5 and 11 sec in travel time x. 

The data presented in Tables 1 through 3 demonstrate that the two algorithms are 
in verv’ good agreement. Note that the accuracy of the RRA algorithm can always be 
improved by reducing the arc length step size. This is demonstrated in Tables 5 and 6. 
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Table 4. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3 



.Method 


>0 (m) 


/^o (Jeg) 


r (km) 




P (deg) 


T (sec) 


s (km) 


ODE 


100.0 


85.0 


10.0 


37.06 


87.244 


6.710387 


10.0183 


RR/\ 


100.0 


85.0 


10.0 


36.63 


87.224 


6.710335 


10.0184 


ODE 


100.0 


95.0 


10.0 


162.94 


92.756 


6.710387 


10.0183 


RRA 


100.0 


95.0 


10.0 


163.37 


92.776 


6.710335 


10.0184 


ODE 


100.0 


135.0 


10.0 


144.9^ 


134.922 


9.448852 


14.1104 


RRA 


100.0 


135.0 


10.0 


144.24 


134.930 


9.449238 


14.1110 



For SSP 4 (parabolic SSP), Table 4 shows that at a horizontal range r = 10 km, with 
the ocean bottom at 200 m, the dilTerences between the two algorithms are as follows: 

• between 0.43 and 0.73 m in depth , 

• between 0.008 and 0.02 deg in angle of arrival ^ , 

• between 52 and 386 n sec in travel time t , and 

• between 0.1 and 0.6 m in path length s. 

As is indicated in Table 4, the differences in travel time are relatively large when 
compared with previous cases. This lower level of accuracy was considered unacceptable 
since, in general, accuracy in the millisecond region is normally required for travel time. 
An attempt to improve accuracy was made by reducing the arc length increment used 
in the RRA algorithm from 2.0 to 1.0 m (see Tables 5 and 6). 

The parabolic SSP which produced Table 4 has the most extreme gradient of all the 
SSPs examined and, thus, the two rays launched at 85° and 95° proceed down range 
without reflection from the surface or bottom, as illustrated in Figures 3 and 4. This 
extreme gradient was considered to have been the probable cause of the discrepancies 
noted above and, therefore, the gradient was reduced and the ocean depth increased to 
permit an examination of a similar profile with a more typical gradient (see Table 7). 
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RECURSIVE RRY RCOUSTICS 




Figure 3. SSP 4 - No Reflections for /5o = 85 ° . 
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RECURSIVE RflY RCOUSTICS 




Figure 4. SSP 4 - No Reflections for Pq= 95° . 
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Table 5. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 2.0 
M ARC LENGTH STEP SIZE 



.Method 


>0 (m) 


/?0 (deg) 


r (km) 


y (m) 


p (deg) 


T (sec) 


s (km) 


ODE 


100.0 


85.0 


5.0 


64.25 


85.597 


3.354779 


5.0088 


RRA 


100.0 


85.0 


5.0 


64.13 


85.581 


3.354757 


5.0088 


ODE 


100.0 


95.0 


5.0 


135.75 


94.403 


3.354779 


5.0088 


RIGA 


100.0 


95.0 


5.0 


135.87 


94.419 


3.354757 


5.0088 


ODE 


100.0 


135.0 


5.0 


77.50 


45.019 


4.724411 


7.0552 


RRA 


100.0 


135.0 


5.0 


77.66 


45.017 


4.724506 


7.0553 



The data in Table 5 for a horizontal range of 5 km and using an arc length step size 
of 2.0 m was generated in order to enable a direct comparison with the arc length step 
size in the RR.A algorithm set to 1.0 m, as shown in Table 6. The reduction in final 
range from 10 to 5 km was necessitated by the structure of the RIGA computer code, 
which required large quantities of computer storage space to be allocated, particularly 
as arc length step size was decreased. This drawback could be rectified by straightfor- 
ward modifications to the simulation program. 

For SSP 4 (parabolic SSP), Table 5 shows that at a horizontal range of r = 5 km, 
the differences between the two algorithms are as follows; 

• between 0.12 and 0.16 m in depth , 

• between 0.002 and 0.016 deg in angle of arrival ^ , 

• between 22 and 95 n sec in travel time t , and 

• 0.1 m in path length ^ 
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Table 6. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 1.0 
M ARC LENGTH STEP SIZE 



Method 


Jo (m) 


Po (deg) 


r (km) 


V (m) 


P (deg) 


T (sec) 


5 (km) 


ODE 


100.0 


85.0 


5.0 


64.25 


85.597 


3.354779 


5.0088 


RRA 


lOO.O 


85.0 


5.0 


64.19 


85.589 


3.354768 


5.0088 


ODE 


100.0 


95.0 


5.0 


135.75 


94.403 


3.354779 


5.0088 


RRA 


lOO.O 


95.0 


5.0 


135.81 


94.411 


3.354768 


5.0088 


ODE 


100.0 


135.0 


5.0 


77.50 


45.019 


4.724411 


7.0552 


RRA 


100.0 


135.0 


5.0 


77.56 


45.018 


4.724458 


7.0553 



For SSP 4 (parabolic SSP), Table 6 shows that at a horizontal range of r = 5 km, 
and with the RRA arc length step size set to 1.0 m, the diflerences betw’een the two al- 
gorithms are as follows: 

• 0.06 m in depth , 

• between 0.00 1 and 0.008 deg in angle of arrival p , 

• between 1 1 and 47 ix sec in travel time t , and 

• 0. 1 m in path length s. 

A significant improvement in the consistency of travel time calculations can be seen 
by comparing Tables 5 and 6. Cutting the arc length increment in half reduced the travel 
time discrepancy by approximately one-half. Further reductions in arc length increment, 
an input parameter to the simulation program, would yield further improvements in 
travel time accuracy. 
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Table 7. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 2.0 
M ARC LENGTH STEP SIZE, BOTTOM DEPTH 1000 M 



Method 


>0 (m) 


Po (deg) 


r (km) 


y (m) 


P (deg) 


T (sec) 


s (km) 


ODE 


500.0 


85.0 


10.0 


177.41 


93.388 


6.705898 


10.0150 


RICA 


500.0 


85.0 


10.0 


777.48 


93.389 


6.705893 


10.0150 


ODE 


500.0 


95.0 


10.0 


222.59 


86.612 


6.705898 


10.0150 


RRA 


500.0 


95.0 


10.0 


222.52 


86.611 


6.705893 


10.0150 


ODE 


500.0 


135.0 


10.0 


545.01 


134.997 


9.448814 


14.1104 


RR.*\ 


500.0 


135.0 


10.0 


545.00 


134.997 


9.448829 


14.1104 



For SSP 4 (parabolic SSP) with vertex c(500) = 1490 m/s, speed of sound at the 
ocean surface c(0) = 1500 m/s, and speed of sound at the ocean bottom c(lOOO) = 1500 
m/s, i.e., an ocean depth of 1000 m; with c{y) = c(500) + a(y — 500)^ w'here 
a = 0.00004, at a horizontal range r = 10 km, the differences between the two algo- 
rithms are as follows: 

• between 0.01 and 0.07 m in depth 

• 0.001 deg in angle of arrival p and, 

• between 5 and 15 sec in travel time t. 

Table 7, when compared with Table 4, both pertaining to an arc length step size of 
2.0 m, demonstrates that a more realistic gradient of 0.02 1/sec enables the RRA algo- 
rithm to handle propagation in a SOFAR channel to a range of 10 km with a much 
improved degree of accuracy. Of course, this could be improved still further by reducing 
the arc length step size. 
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B. RRA THREE-DIMENSIONAL CAPABILITIES 



Figure 5 demonstrates the ability of the RRA algorithm to handle launch angles 
other than 0,, = 90° . This figure u'as produced using the parabolic SSP 4 with 

c(j’) = cOo) + «0 - Vo)". 0 < >>< 200m, (3.1) 

where j’o = 100 m, c(vo) = 1490.0 m/s, a = 0.001 ^ , and </)o = 85° . Note that the 

ray stays in the plane created by this angle, as illustrated by the projection of the ray 
onto the X Z plane as a straight line. 

The RRA algorithm also has the capability to divide the ocean into separate zones 
as a function of z. Each zone can have its own SSP as a function of both cross range jc 
and depth y. For computer simulation purposes, the speed of sound was defined as fol- 
lows: 

c{x, >-) = c(ko) + aiy-y^f + , 0 < y < 200 m , 

where yo = 100 m, c(y’o) = 1490.0 m/s, a = 0.001 ^ and 

r 0.0 1/s, z < 5 km 
j-0.05 1/s, 5 < z < 10 km 

The particular gradient in the x direction was selected solely for computer simulation 
purposes to illustrate the capability of the RRA algorithm to acconvmodate sound speed 
as a function of both jr and y. Figure 6 illustrates the ray trace provided by the RRA 
algorithm for the SSP given by Equations 3.2 and 3.3, showing the ray staying in the 
Y Z plane until reaching 5 km, where the gradient in the jc direction comes into effect, 
resulting in out-of-plane propagation. 



(3.2) 



(3.3) 
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Figure 6. Ray Trace Based on a SSP as a Function of x, y and z • 
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IV. REPRESENTATION OF SOUND-SPEED DATA 



A. COMPARATIVE ANALYSIS - SPATIAL FOURIER SERIES VERSUS AKIMA 
CUBIC SPLINE REPRESENTATIONS OF SOUND-SPEED DATA 

Five diOercnt exact mathematical functions were used to generate II sound-speed 
data points evenly spaced in depth in a 200 m deep ocean. Akima cubic spline and 
spatial Fourier series representations of the SSPs were produced from these data points. 
Percent error calculations of their sound-speed fits as well as their first and second-order 
derivative fits were performed with respect to the known exact values. If the error be- 
tween the exact and modeled values was less than 10~‘ , percent error was set to zero. 

The exact functions examined were as follows: 

• Constant speed of sound ( c = 1500 m's ) 

• Linear SSP with a positive gradient ( c(j') = 1500.0 + O.OITj; m/s ) 

• Linear SSP with a negative gradient ( c(j-) = 1500.0 — O.OlTy m/s ) 

/tv 

• Half-period sine wave SSP ( c(/v) = 1500.0 — 25.0 sin( ) m/s ) 

• Parabolic SSP ( dj’) = 1450.0 + 0.005Cv - 100.0)' m/s ) 

Tables 8 and 9 with Figures 7, 8 and 9 illustrate the dilTiculties that the Fourier 
representation encounters in dealing with data derived from an exact constant speed of 
sound. No noticeable difference in the sound-speed It itself is observed for either the 
Akima cubic spline or the Fourier series representation, but the latter generates some 
small errors in both the first and second-order derivatives. 

Tables 10 and 11 along with Figures 10, 11 and 12 compare the Akima cubic spline 
and Fourier series representations of sound-speed data generated from a linear SSP with 
a positive gradient. In Figure 10 and Table II it can be seen that the Fourier represen- , 
ration does not precisely match the sound-speed function, although Figure 10 and Table 
10 show that the Akima cubic spline does. Problems similar to those noted for a con- 
stant speed of sound can be observed by examining Figures II and 12 and Table II. 
Discrepancies in the derivatives are now much larger. 

Tables 12 and 13 along with Figures 13, 14 and 15 show analogous results for a 
linear SSP with a negative gradient. Once again, neither the sound-speed fit nor the 
derivatives are precisely modeled by the Fourier series representation. 
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Table 8. AKIMA CUBIC SPLINE RESULTS - CON- 
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NUMBER OF ORIGINAL DATA POINTS 



Table 9 



FOURIER SERIES RESULTS - CONSTANT 
SPEED OF SOUND 
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Figure 8, First-Order Derivative Fit - Constant Speed of Sound. 
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Figure 9. 



Second-Order Derivative Fit - Constant Speed of Sound. 
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Table 10, AKIMA CUBIC SPLINE RESULTS - LINEAR 
SSP WllH A POSITIVE GRADIENT 
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Table 11. FOURIER SERIES RESULTS - LINEAR SSP 
WITH A POSITIVE GRADIENT 
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Figure 11. First-Order Derivative Fit - Linear SSP with a Positive Gradient. 
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TECHNIQUES FOR NUMERICRL INTERPOLnTION OF C(Y) 




Figure 12. Second Order Derivative Fit - Linear SSP with a Positive Gradient. 
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Table 12. AKIMA CUBIC SPLINE RESULTS - LINEAR 
SSP WITH A NEGATIVE GRADIENT 
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Table 13. FOURIER SERIES RESULTS - LINEAR SSP 
WITH A NEGATIVE GRADIENT 
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Figure 13. Sound Speed Fit - Linear SSP with a Negative Gradient. 
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LINERR SSP WITH fi NOTTIVE BRfiOIENT 

ORIOINRL NUMBER OF DRIB POINTSill NUTIBER OT SYNTHETIC DRTfl POINTS BEFORE fM) RFTERiO 

WITHOUT UlOOS SnOOTHlNQ 




Figure 1*4. 



First-Order Derivative Fit - Linear SSP with a Negative Gradient. 
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Figure 15. Second-Order Derivative Fit - Linear SSP with a Negative Gradient. 
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ONIHIOOUS SOZONtn inOHllM 



The above three exact sound-speed functions were used to compare the Fourier se- 
ries representation with the Akima cubic spline since the Fourier fit can be expected to 
be at its worst since we are tr\'ing to form linear functions from trigonometric ones. The 
next two functions examined, however, are non-linear; the half-period sine wave and the 
parabola. 

Tables 14 and 15 along with Figures 16, 17 and 18 compare the Akima cubic spline 
and Fourier series representations of sound-speed data generated from a half-period sine 
wave. Both techniques encounter some difficulty here, with the Fourier representation 
generally behaving better near the midpoint of the data than the Akima cubic spline fit 
for sound speed and both derivatives. However, near the surface and bottom, the errors 
in the Fourier representation become much larger than those of the Akima cubic spline 
fit. 

Tables 16 and 17 along with Figures 19, 20 and 21 compare the Akima cubic spline 
and Fourier series representations of sound-speed data generated from a parabola. 
Surprisingly, the Akima cubic spline representation matches sound speed and derivatives 
ver>' well. The Fourier representation behaves in a manner siinilar to its handling of the 
half- period sine wave. It produces small errors near the midpoint of the data and large 
discrepancies at either end for sound speed as well as in the first and second-order de- 
rivatives. 

The disappointing performance of the Fourier series representation of sound speed 
and of the first and second-order derivatives of sound speed led to various attempts to 
salvage it. First, Lanczos smoothing as discussed in Chapter II, Section C, w'as imple- 
mented. This ameliorated the oscillations of the Fourier representation as shown in 
Tables 18 and 19 and Figures 22 through 27 for the half-period sine wave and parabolic 
SSPs. It can be seen that Lanczos smoothing results in a smoother fit for sound speed 
and for the derivatives to a point closer to the end points of the data but, nevertheless, 
leaves significant discrepancies near the end points. 
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TECHNIQUE! AKIMA CUBIC SPLINE FIT 



Table 14. AKIMA CUBIC SPLINE RESULTS 
HALF-PERIOD SINE WAVE SSP 
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Table 15. FOURIER SERIES RESULTS - HALF-PERIOD 
SINE WAVE SSP 
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TECHNIQUES FOR NUMERICRL INTERPOLRTION OE C(Y) 




Figure 16. Sound Speed Fit - Half-Period Sine Wave SSP. 
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TECHNIQUES FOR NUMERICRL INTERPOLRTION OF C(Y) 




Figure 18 . Second-Order Derivative Fit - Half-Period Sine Wave SSP. 
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Table 17. FOURIER SERIES RESULTS - PARABOLIC 
SSP 
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TECHNIQUES FOR NUhRRICRL INTERPOLRTION OF C(Y) 




Figure 19. Sound Speed Fit - Parabolic SSP. 
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6NIH1QOUS SOZOHn inOHlIH 




Figure 20. First-Order Derivative Fit - Parabolic SSP. 
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Figure 21. 



Second-Order Derivative Fit - Parabolic SSP. 
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Table 



18. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING - HALF-PERIOD SINE WAVE 
SSP 
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NUMBER OF ORIGINAL DATA POINTS = 11 

TECHNIQUE: SPATIAL FOURIER SERIES FIT WITH LANCZOS SMOOTHING 



Table 19. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING - PARABOLIC SSP , 
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Figure 22. Sound Speed Fit with Lauczos Smoothing - Half-Period Sine Wave 
SSP. 
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'igure 23. 



First-Oider Derivative Fit with Lanczos Smoothing 



Half- Period Sine 



Wave SSP. 
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Figure 24. Second-Order Derivative Fit witli Lanczos Smoothing - Half-Period Sine 



Wave SSP. 



55 



TECHNIQUES FOR NUNERICni INTERPOLATION OF C(Y) 




Figure 25. Sound Speed Fit with Lanczos Smoothing - Parabolic SSP. 
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ORISINBL NUMBER OT DRTfl POINrSill NUMBER DP SYNTHETIC DniR POINTS BEFORE AND fFIERiO 

WITH LflNCZOS SMOOTHINO 
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Figure 26. First-Order Derivative Fit with Lanczos Smoothing - Parabolic SSP. 
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TECHNIQUES FOR NUMERICRL INTCRPOLRTION OF C(Y) 
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Figure 27. Second-Order Derivative Fit uilli Lanczos Smoothing - Parabolic SSP. 
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ORICINRL NUMBER OT DRTfl POINTSill NUMBER OP SYNTHETIC DftTR POINTS BEFORE RND flPTERiO 

WITH LRNCZOS SM00IHIN3 



The difficulty in achieving a good fit at or near the end points of the data using a 
spatial Fourier series representation is a result of the instantaneous commencement and 
cessation of sound-speed data at the ocean surface and bottom. Therefore, an attempt 
to overcome this difficulty was made by adding synthetic data points above the ocean 
surface and below the ocean bottom so as to enable the Fourier series to "settle dovvn" 
before the real data begins. 

Synthetic data was created by first examining the three data points nearest the ocean 
surface and bottom. Similar procedures were used to generate synthetic data both above 
the surface and below the bottom. Therefore, only the surface case will be exainined. 

Let 



m 



10 



cfvi) - c(0) 
Ti 



where c(0) = c(y = 0) , and 



cCv'2) - 

yi >1 



where is depth in meters at the ith data point. If we let 

A/n = /« 2 i — 



and if 



|A/?i| = \ni 21 - m,oI <t 



for some arbitrarily small e, then linear extrapolation is used, i.e., 

cO-i) = c(0) + i= \,2,-,SYNPTS 

where 



y~i = -yi 

and SYNPTS is the number of synthetic data points desired. 



(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
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If the test at Equation 4.4 fails, and if 

■ "^lol < l"»2i I - (4.7) 

then mirror imaging about the ocean surface is appropriate, i.e., 

c(y-i) = cO;), / = l,2,...,SYNPTS . (4.8) 

If both of the tests at Equations 4.4 and 4.7 fail, then neither linear extrapolation 
nor mirror imaging is appropriate. Under these circumstances it is appropriate to re- 
verse the rate of change of the slope and let 

cOJ = c(0) -f (m,o - iAm)y_i, i = \,2,-^,SYNPTS . (4.9) 

This methodology of extending the real data was implemented using five synthetic 
data points both above the surface and below the bottom. The results are shown in 
Tables 20 through 24 and Figures 28 through 42. 

A substantial improvement in the accuracy of the Fourier series representation can 
be seen with respect to all the functions examined. Nevertheless, some difficulty is still 
evident near the ends of the original data. A final attempt to resolve this discrepancy 
was made by using the Akima cubic spline representation of sound speed as input to the 
Fourier algorithm, and sampling the SSP generated by the Akima cubic spline fit at some 
arbitrarily high frequency. This was attempted using a sampling rate ten times greater 
than available from the original data. The synthetic data as described above was not fit 
using the Akima cubic spline method. Results for the half-period sine wave and 
parabolic SSPs are presented in Tables 25 and 26 and Figures 43 to 48. 

These results are once again an improvement on previous work. Nonetheless, some 
inaccuracies remain. Unfortunately, time constraints precluded further, more sophisti- 
cated attempts to further improve the accuracy of the Fourier representation. 
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Table 20. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING AND SYNTHETIC DATA - 
CONSTANT SPEED OF SOUND 
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Table 21. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING AND SYNTHETIC DATA - LIN- 
EAR SSP WITH POSn IVE GRADIENT 
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Table 22. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING AND SYNTHETIC DATA - LIN- 
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Table 23. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING AND SYNTHETIC DATA - 
HALF-PERIOD SINE WAVE SSP 
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Table 24. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING AND SYNTHETIC DATA - 
PARABOLIC SSP 
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TECHNIQUES FOR NUMERICRL INTERPOLRTIQN OF C(Y) 




Figure 28. Sound-Speed Fit - Constant Speed of Sound with Synthetic Data. 
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ORIOIKRL NUilDER OT DfTTR POINISill NUMBER OF SYNTHETIC OflTfi POINTS BEFORE RNO fiFTERiS 

WITH LRNC2DS SMOOTHINQ 
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Figure 29. First-Order Derivative Fit 
Data. 



Constant Speed of Sound nith Syntlietic 



67 



ORIOIMHL MIMDER OT OfTTfl POINTSilI NUMEEIR OT SWnOIC ORTfl POINTS BEFORE fWD fTTERiS 

WITH IBNCZDS SMOOTHIMS 
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Figure 30. Second-Order Derivative Fit - Constant Speed of Sound with Synthetic 



Data. 
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TECHNIQUES FOR NUHERICRL INTCRPOLFITION Or C(Y) 




Figure 31. Sound-Speed Fit - Linear SSP with a Positive Gradient and Synthetic 
Data. 
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ORIQINRL NUMBER OT Dftm POINTSill NUMBER OT SYNTHETIC DBTR POINTS BEFORE RNO RHERiS 

WITH LniiCZOS SMOOTHna 
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Figure 32. 



First-Order Derivative Fit 



Linear SSP witii a Positive Gradient and 



Synthetic Data. 
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Synthetic Data. 
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Figure 34. Sound-Speed Fit - Linear SSP with a Negative Gradient and Synthetic 
Data. 
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Figure 35. First-Order Derivative Fit - Linear SSP nitli a Negative Gradient and 
Synthetic Data. 
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Figure 36. Second-Order Derivative Fit - I inear SSP with a Negative Gradient and 



Synthetic Data. 
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Figure 38. First-Order Derivative Fit - Half-Period Sine Wave SSP with Synthetic 



Data. 
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TECHNIQUES FOR NUhERICRL INTERPOLRTION OF C(Y) 
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Figure 39. Second-Order Derivative Fit - Half-Period Sine Wave SSP with Syn- 



thetic Data. 
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0RI81NRL KUMBE3? Of ORTfl POlNTSill NIH5CR OT SYWHOIC DflTfl POINTS BtTORE AND ATTERiS 

WITH LnNCZOS SM00THIII3 



TECHNIQUES FOR NUhERICRL INTERPOLRTION OF C(Y) 




Figure 40. Sound-Speed Fit - Parabolic SSP with Synthetic Data. 
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ORIOINRL NUMBER OT DRTfl POINTSill NUMBER OF SOTHETIC OflTfl POINTS BEFORE fWD RFTERiS 

WITH LRNCZOS SMOOTHINO 
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Figure 41. First-Order Derivative Fit - Parabolic SSP with Synthetic Data. 
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ORIGINRL NUH0E3^ OF DRTfl POINTSill KUMBER OF SYNTHniC DflTR POINTS BEFORE RND FTTERiS 
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Figure 42. Second-Order Derivative Fit - Parabolic SSP with Synthetic Data. 
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Table 25. FOURIER SERIES RESULTS WITH LANCZOS 
SMOOTHING, SYNTHETIC DATA, AKIMA 
CUBIC SPLINE INPUT TO FOURIER 
METHOD - HALF-PERIOD SINE WAVE SSP 
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Table 26 . FOURIER SERIES RESULTS WITH LANCZOS SMOOTHING, 
SYNTHETIC DATA, AKIMA CUBIC SPLINE INPUT TO FOURIER 
METHOD - PARABOLIC SSP 
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TECHNIQUES FOR NUhERICRL INTERPOLATION OF C(Y) 




Figure 43. Sound-Speed Fit - Half-Period Sine Wave SSP, Akiina Cubic Spline 
Input to Fourier Method. 
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ORIOINnL MJMDES GF DRTR POINTSiII NUMBCR OT ST^mOIC DflTfl POINTS BCTORE fiNO fiFTERiS 

WITH LflNCZOS SMOOTHINS SPflTinL FOURIER SERIES FIT OBTfllNEO BY Sfm-ING mm CUBIC SPLINE FIT WITH SflMPLINQ RBTE 10 




Figure 44. 



First-Order Derivative Fit - Half-Period Sine Wave SSP, Akima Cubic 
Spline Input to Fourier Method. 
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"igure 45. 




Second-Order Derivative Fit - Half-Period Sine Wave SSP, Akima Cu- 



bic Spline Input to Fourier Method, 
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TECHNIQUES FOR NUMERICRL INTERPOLRTION OF C(Y) 




Figure 47. First-Order Derivative Fit - Parabolic SSP, Akitna Cubic Spline Input 
to Fourier Method. 
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Figure 48. Second-Order Derivative Fit - Parabolic SSP, AkJma Cubic Spline Input 



to Fourier Method. 
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B. RAY TRACE RESULTS USING SPATIAL FOURIER SERIES 
REPRESENTATION OF SOUND-SPEED DATA WITH THE RRA ALGORITHM 
1. Sound Speed Including Deterministic Component Only 

The RR/\ algorithm using the simple Fourier series sound-speed representation 
including Lanczos smoothing but without synthetic data or Akinia cubic spline input to 
the Fourier method was tested using the same source locations, SSPs, and rays as in 

J 

Chapter III, Section A. That is, the SSPs were once again defined as follows: 

• SSP 1: Constant speed of sound (c = ISOOmjs) 

• SSP 2: Linear SSP with a positive gradient ( c{y) = c(0) -f gy m/s, where the 
speed of sound at the surface c(0) = c(y = 0) = 1500 m/s and g = -t-0.017 1/s ) 

• SSP 3; Linear SSP with a negative gradient ( = c(0) + gy m/s, where the 
speed of sound at the surface c(0) = c(y = 0) = 1500 m/s and g = —0.017 1/s ) > 

• SSP 4: Parabolic SSP with vertex at y^ where c(j’o) = 1490 m/s and the speed of 
sound at both the ocean surface and bottom is equal to 1500 m,s 

(c(y) = c(yo) + a (y — 100)^ m/s with a = 0.001 for an ocean depth of 200 

m ). 

The sound source spatial location was established at .r,, = 0 m and z,, = 0 m. A 
source depth ofy^ = 100 m was established for an ocean depth of 200 m, unless otherwise 
noted. Launch angles /?o were 45®, 85® and 135® for SSPs 1, 2 and 3; and 85®, 95®, and 
135® for SSP 4, which again was designed to simulate a SOFAR channel. 

The ODE algorithm using the Akima cubic spline fit to sound-speed data re- 
mained as the benchmark for comparison purposes, as detailed in Chapter III, Section 
A. For ease of reference, the results of the ODE algorithm are reproduced in the tables 
of this sub-section. The ODE's range step size was left at 2.0 m and the arc length step 
size in the RRA algorithm was left at 2.0 m as well, except where otherwise noted. 

The results of the various computer simulation runs are presented in Tables 27 
to 33 below. 
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Table 27. RAY ACOUSTICS RESULTS FOR SSP 1, RAYS 1 THROUbH 3, 
FOURIER SOUND-SPEED FIT 



Method 


A (m) 


Po (deg) 


r (km) 


V (m) 


P (deg) 


T (sec) 


s (km) 


ODE 


100.0 


45.0 


10.0 


100.00 


45.000 


9.428090 


14.1421 


RR.A 


100.0 


45.0 


10.0 


99.94 


45.000 


9.428088 


14.1421 


ODE 


100.0 


85.0 


10.0 


174.89 


85.000 


6.692132 


10.0382 


RR.A 


100.0 


85.0 


10.0 


174.88 


85.000 


6.692132 


10.0382 


ODE 


100.0 


135.0 


10.0 


100.00 


135.000 


9.428090 


14.1421 


RRA 


100.0 


135.0 


10.0 


100.06 


135.000 


9.428088 


14.1421 



For SSP 1 (constant speed of sound), Table 27 shows that at a horizontal range 
of r= 10 km, the differences between the benchmark ODE algorithm and the RRA al- 
gorithm using the Fourier series sound-speed fit are as follows: 

• between 0.01 and 0.06 m in depth and 

• 2 yu sec in travel time t. 

Note the significant degradation of the RRA results relative to those obtained 
using the Akima cubic spline fit as in Table 1 of Chapter III, Section A. This compar- 
ative degradation of consistency between the ODE and RRA algorithms will be evident 
throughout this subsection. 
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Table 28. RAY ACOUSTICS RESULTS FOR SSP 2, RAYS 1 THROUGH 3, 
FOURIER SOUND-SPEED FIT 



.Method 


do (m) 


Po (deg) 


r (km) 


d (m) 


P (deg) 


T (sec) 


s (km) 


ODE 


100.0 


45.0 


10.0 


99.99 


45.000 


9.417413 


14.1421 


RRA 


100.0 


45.0 


10.0 


99.97 


45.000 


9.417429 


14.1421 


ODE 


100.0 


85.0 


10.0 


162.5C 


85.487 


6.683492 


10.0374 


RRA 


100.0 


85.0 


10.0 


163.31 


85.482 


6.683550 


10.0374 


ODE 


100.0 


135.0 


10.0 


100.01 


135.000 


9.417413 


14.1421 


RRA 


100.0 


85.0 


10.0 


100.03 


135.000 


9.417430 


14.1422 



For SSP 2 (linear SSP with a positive gradient), Table 28 shows that at a hori- 
zontal range r = 10 km, the difTerences between the benchmark ODE algorithm and the 
RRA algorithm using the Fourier series sound-speed fit are as follows; 

• between 0.02 and 0.19 m in depth 

• 0.005 deg in angle of arrival P, 

• between 16 and 58 fx sec in travel time t, and 

• 0.1 m in path length 5. 
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Table 29. RAY ACOUSTICS RESULIS FOR SSP 3, RAYS 1 THROUGH 3, 
FOURIER SOUND-SPEED FIT 



Method 


To (m) 


Po (deg) 


r (km) 


y (nt) 


P (deg) 


T (sec) 


s (km) 


ODE 


100.0 


45.0 


10.0 


99.99 


45.000 


9.438784 


14.1421 


RRA 


100.0 


45.0 


10.0 


99.96 


45.000 


9.438799 


14.1422 


ODE 


100.0 


85.0 


10.0 


168.89 


84.512 


6.699404 


10.0380 


RRA 


100.0 


85.0 


10.0 


169.76 


84.492 


6.699471 


10.0380 


ODE 


100.0 


135.0 


10.0 


100.01 


135.000 


9.438784 


14.1421 


RIU\ 


100.0 


85.0 


10.0 


100.04 


135.000 


9.438798 


14.1421 



For SSP 3 (linear SSP with a negative gradient), Table 29 shows that at a hori- 
zontal range r = 10 km, the difference between the benchmark ODE algorithm and the 
RRi\ algorithm using the Fourier sound-speed fit are as follows: 

• between 0.03 and 0.13 m in final depth j;, 

• 0.020 deg in angle of arrival /?, 

• between 14 and 67 n sec in travel time t, and 

• 0.1 m in path length s. 

The results contained in the above tables, while poorer than those detailed in 
Chapter III, Section A, nevertheless produce travel times which are consistent in the 
millisecond region. However, this is not the case for the SSPs examined next. 
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Table 30. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 
FOURIER SOUND-SPEED FIT 



Method 


;-o (m) 


Po (deg) 


r (km) 


.V (m) 


P (deg) 


T (sec) 


s (km) 


ODE 


100.0 


45.0 


10.0 


37.06 


87.244 


6.710387 


10.0183 


RRA 


100.0 


45.0 


10.0 


42.86 


86.674 


6.709733 


10.0185 


ODE 


100.0 


85.0 


10.0 


162.94 


92.756 


6.710387 


10.0183 


RRy\ 


100.0 


85.0 


10.0 


157.14. 


93.326 


6.709733 


10.0185 


ODE 


100.0 


135.0 


10.0 


144.97 


134.922 


9.448852 


14.1104 


RRA 


100.0 


85.0 


10.0 


143.98 


134.930 


9.448515 


14.1112 



For SSP 4 (parabolic SSP), Table 30 shows that at a horizontal range of r = 10 
km, the differences between the benchmark ODE algorithm and the RRA algorithm 
using the Fourier series sound-speed fit are as follows: 

• between 0.99 and 5.8 m in final depths, 

• between 0.008 and 0.570 deg in angle of arrival P, 

• between 337 and 654 n sec in travel time t, and 

• 0.2 m in path length 5. 

As in Chapter III, Section A, SSP 4 shows significantly less consistency be- 
tween the two algorithms than the previous SSPs examined. Again, two approaches 
were taken to attempt to reduce this discrepancy. First, the arc length step size was re- 
duced from 2.0 to 1.0 m (see Tables 12 and 13) and second, the steepness of the gradient 
was reduced from 0.1 to 0.02 1/sec (see Table 14). 
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Table 31. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 2.0 
M ARC LENGTH STEP SIZE, FOURIER SOUND-SPEED FIT 



Method 


Vo (m) 


k (deg) 


r (km) 


V (m) 


/? (deg) 


T (sec) 


5 (km) 


ODE 


100.0 


45.0 


5.0 


64.25 


85.597 


3.354779 


5.0088 


RIU\ 


100.0 


45.0 


5.0 


68.79 


85.435 


3.354591 


5.0090 


ODE 


100.0 


85.0 


5.0 


135.75 


94.403 


3.354779 


5.0088 


RRA 


100.0 


85.0 


5.0 


131.21 


94.565 


3.354591 


5.0090 


ODE 


100.0 


135.0 


5.0 


77.50 


45.019 


4.724411 


7.0552 


RR.A 


100.0 


85.0 


5.0 


77.82 


45.017 


4.724162 


7.0554 



For SSP 4 (parabolic sound speed), Table 31 shows that at a horizontal range 
of r = 5 km, using an arc length step size of 2.0 m, the differences between the 
benchmark ODE algorithm and the R1U\ algorithm using the Fourier series sound-speed 
fit are as follows; 

• between 0.32 and 4.54 m in final depth 

• between 0.002 and 0.162 deg in angle of arrival /J, 

• between 188 and 249 ^JL sec in travel time t, and 

• 0.2 m in path length s. 

The data in Table 31 above was generated in order to enable a direct compar- 
ison with the arc length step size set to 1.0 m, as shown in Table 32 below. The re- 
duction in final range from 10 to 5 km was necessitated by the structure of the RRA 
computer code, which required large quantities of computer storage space to be allo- 
cated, particularly as arc length step size was reduced. 
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Table 32, RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 1.0 
M ARC LENGTH STEP SIZE 



Method 


>0 (m) 


Po (deg) 


r (km) 


>’{m) 


/? (deg) 


T (sec) 


5 (km) 


ODE 


100.0 


45.0 


5.0 


64.25 


85.597 


3.354779 


5.0088 


RRA 


100.0 


45.0 


5.0 


68.86 


85.442 


3.354604 


5.0090 


ODE 


100.0 


85.0 


5.0 


135.75 


94.403 


3.354779 


5.0088 


RRi\ 


100.0 


85.0 


5.0 


131.14 


94.558 


3.354604 


5.0090 


ODE 


100.0 


135.0 


5.0 


77.50 


45.019 


4.724411 


7.0552 


RRj\ 


100.0 


85.0 


5.0 


77.74 


45.019 


4.724126 


7.0554 



For SSP 4 (parabolic SSP), Table 32 shows that at a horizontal range of r = 5 
km, and with the RRA arc length step size set to 1.0 m, the differences between the 
benchmark ODE algorithm and the RRA algorithm using the Fourier series sound-speed 
fit are as follows: 

• between 0.24 and 4.61 m in final depth j;, 

• 0.155 deg in angle of arrival /?, 

• between 175 and 285 ix sec in travel time t, and 

• 0.2 m in path length s. 
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Table 33. RAY ACOUSTICS RESULTS FOR SSP 4, RAYS 1 THROUGH 3, 2.0 
M ARC LENGTH STEP SIZE, BOTTOM DEPTH 1000 M, FOURIER 
SOUND SPEED FIT 



.Method 


>'o (m) 


/Jo (deg) 


r (km) 


>'(ni) 


/? (deg) 


T (sec) 


5 (km) 


ODE 


100.0 


45.0 


10.0 


777.41 


93.388 


6.705898 


10.0150 


RR.A 


100.0 


45.0 


10.0 


768.46 


93.522 


6.705420 


10.0152 


ODE 


100.0 


85.0 


10.0 


222.59 


86.612 


6.705898 


10.0150 


RR/\ 


100.0 


85.0 


10.0 


231.54 


86.478 


6.705420 


10.0152 


ODE 


100.0 


135.0 


10.0 


545.01 


134,997 


9.448829 


14.1104 


RRA 


100.0 


85.0 


10.0 


544.58 


134.997 


9.448180 


14.1107 



For SSP 4 (parabolic speed of sound), with vertex c(500) = 1490 m/s and speed 
of sound at the ocean surface c(0) = 1500 m/s and speed of sound at the ocean bottom 
c(lOOO) = 1500 m/s, i.e., an ocean depth of 1000 m; at a horizontal range of r = 10 km, 
the differences between the benchmark ODE algorithm and the RRA algorithm using 
the Fourier series sound-speed fit are as follows: 

• between 0.43 and 8.95 m in final depth y, 

• 0.134 deg in angle of arrival /?, 

• between 478 and 649 /u sec in travel time t, and 

• between 0.2 and 0.3 m in path length s. 

Table 33 shows no significant improvement over Table 30 for SSP 4 with the 
original extreme gradient. The errors associated with the Fourier sound speed repre- 
sentation have resulted in this lack of improvement, whereas substantial improvement 
was observed between these tw'’o cases for the Akima cubic spline fit with the RRA al- 
gorithm (see Tables 4 and 7 in Chapter III, Section A). 
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2. Sound Speed Including Deterministic and Random Components 

In order to investigate the efTects of randomness in the SSP, a modi^l of the 
power spectral density of the random component of the speed of sound CrO--) had to be 
proposed. The most straightforward model of band-limited white noise was used, as 
given by 




5o, \fy\ < Biv 
0, \fy\ > BIV 



(4.10) 



where BW is the bandwidth, and So is the constant magnitude of the power spectrum at 
all spatial frequencies . 

As a consequence of the assumption of zero mean white noise, the variance of 
c«(j.) is given by 

a]^=2S,BW (4.11) 



and, as a result. 






2 




2BW ■ 



(4.10) 



A value for the variance of C;,(y) in shallow water of 2.5 x 10"’’ was obtained 
from Brekhovskikh and Lysanov [Ref. 8: p. 208 ]. The simulation program was run with 
randomness incorporated as above and for SSP 1 (constant speed of sound). When the 
results were compared with those from the RRA algorithm previously run using the 
Fourier series representation of sound-speed data, differences of a maximum of 0.1 m in 
final depth y, 0.002 deg in angle of arrival /?, and 2 /i sec in travel time t were observed. 
The RRA algorithm was then run with SSP 4 (parabolic SSP in a 200 m deep ocean) 
incorporating randomness and differences of a maximum of 0.22 m in final depth y , 
0.018 deg in angle of arrival p, and 16 /r sec in travel time t were observed. 

The variance was then increased to determine at what magnitude it caused the 
randomness to become significant. Further simulation runs revealed that for SSP 1, 
differences in travel times of greater than 10 n sec needed al^ to be at least of the order 
of 10'^ , four orders of magnitude greater than the estimate given in Brekhovskikh and 
Lysanov [Ref 8: p. 208 [. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The recursive ray acoustics (RFU\) algorithm was shown to be an accurate ray 
acoustics algorithm when used in conjunction with an Akima cubic spline fit to sound- 
speed data. The RRA algorithm can be made as accurate as desired by reducing the arc 
length step size. .Although the ODE algorithm computes amplitudes which the RICA 
algorithm does not, the RRA algorithm was observed to be significantly faster in terms 
of computer execution time. 

The ability of the RICA algorithm to handle sound speed as a function of all three 
spatial coordinates enables its application to real world problems wherein a horizontal 
gradient of sound speed exists, either in the cross-range or down-range direction. This 
could conceivably be the case for either long range ray propagation or for propagation 
in estuaries or across current boundaries where significant horizontal gradients can exist. 

.A spatial Fourier series (SFS) representation of sound-speed data, while attractive 
for reasons of its spectral content and the ease of incorporation of randomness at spe- 
cific harmonics, encountered dilTiculty with regard to providing accurate fits to sound- 
speed data. Several techniques were used to attempt to rectify this, including the use 
of Lanezos smoothing , addition of synthetic data above the ocean surface and below 
the ocean bottom, and sampling of the Akima cubic spline fit to provide data points for 
input to the SFS method. These techniques resulted in some substantial success in the 
reduction of the inaccuracies of the SFS method. However, although the SFS method 
eventually provided accurate fits to the sound-speed data, it still provided inaccurate first 
and second-order derivatives when compared against the results obtained from the 
Akima cubic spline method. 

A comparison was made between the ODE algorithm using an Akima cubic spline 
sound-speed representation and the RRA algorithm using a SFS sound-speed represen- 
tation, with the SFS representation in its most basic form, that is, with none of the 
techniques referred to in the above paragraph incorporated. The RRA/SFS combina- 
tion was found to be wanting in accuracy. 

As was mentioned previously, the RRA algorithm using the Akima cubic spline fit 
to sound-speed data was proven fast and accurate, and can be made as accurate as de- 
sired by reducing the arc length step size. However, investigation is warranted to de- 
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lermine if the RRA algorithm can be generalized to include amplitude calculations, 
especially at focal points and, hence, caustics. 

Throughout the simulations undertaken in this thesis, a zero-slope, flat ocean sur- 
face and bottom was assumed. In order to apply the RRA algorithm to realistic ocean 
boundaries, the ability to accommodate arbitrarily shaped surfaces should be included. 
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